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ABSTRACT 

We combine an analytic model for anisotropic outflows and galaxy formation with numerical sim- 
ulations of large-scale structure and halo formation to study the impact of galactic outflows on the 
evolution of the Intergalactic medium. Using this algorithm, we have simulated the evolution of a 
comoving volume of size (15Mpc) 3 in the ACDM universe. Using an N-body simulation starting 
at redshift z = 24, we follow the formation of 20 000 — 60 000 galaxies and simulate the galactic 
outflows produced by these galaxies, for five outflow opening angles, a — 60°, 90°, 120°, 150°, and 
180° (isotropic outflows). Anisotropic outflows follow the path of least resistance and thus travel 
preferentially into low-density regions, away from cosmological structures (filaments and pancakes) 
where galaxies form. These anisotropic outflows are less likely to overlap with one another than 
isotropic ones. They are also less likely to hit pre-galactic collapsing halos and strip them of their 
gas, preventing a galaxy from forming. Going from 180° to 60°, the number of galaxies that actually 
form doubles, producing twice as many outflows, and these outflows overlap to a lesser extent. As 
a result, the metal volume filling factor of the IGM goes from 8% for isotropic outflows up to 28% 
for anisotropic ones. High-density regions are more efficiently enriched than low density ones (~ 80% 
compared to ~ 20% by volume), even though most enriched regions are low-densities. Increasing the 
anisotropy of outflows increases the extent of enrichment at all densities, low and high. This is in part 
because anisotropic outflows are more numerous. When this effect is factored-out, we find that the 
probability a galaxy will enrich systems at densities up to 10 p is higher for increasingly anisotropic 
outflows. This is interpreted as an effect of the dynamical evolution of the IGM. Anisotropic outflows 
expand preferentially into underdense gas, but that gas can later accrete onto overdense structures. 
The inclusion of photoionization suppression of low-mass galaxy formation reduces the degree of late 
galaxy formation and preferentially suppresses galaxy formation in low-density regions. The result is 
a decline in the physical extent of galactic outflows after z = 3 as accretion overwhelms the expansion 
of new outflows and reduces feedback in underdense regions. 

Subject headings: cosmology — galaxies: formation — Intergalactic medium — methods: numerical 



II. 



1. INTRODUCTION 

The evolution of the Intergalactic medium (IGM) can 
be significantly affected by feedback effects from galaxies. 
Supernovae and active galactic nuclei (AGN) can deposit 
large amounts of energy into the surrounding interstellar 
gas, accelerating it to large velocities. If the gas becomes 
unbound, a galactic wind, or outflow, will result, with 
important consequences for the evolution of the IGM. 
Galactic outflows deposit energy, momentum, and metal- 
enriched gas into the IGM. This can affect the subsequent 
formation of other galaxies, by striping the gas from col- 
lapsing halos, reheating and possibly ionizing the IGM, 
and modifying its cooling rate through metal enrich- 
ment. These m etals have been observed via the Lyman- 
a forest (e.g. iMever fe Yorkl [l987t ISchave et all 120031: 



Pieri & Hachnelt 
2008t iPieri et al 
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2004 IPieri et all l2009t lAguirre et al.l 
2010). Feedback from galactic out- 



flows have been invoked to explain several observations 
at the galaxy scale, such as the high mass-to-light ratio 
of dwarf galaxies, the overcooling problem, and the size 
of the galactic discs (the angular momentum problem). 
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At larger scales, feedback is needed to explain the ob- 
served metallicity and entropy content of the IGM, and 
the scaling relations for clusters. 

SNe-driven outflows and AGN-driven outflows are 
complementary. The number of SNe increases with 
galactic mass, but the depth of the potential well that 
the gas must climb to escape also increases with galactic 
mass. Hence, outflows produced by high-mass galaxies 
are not much larger than outflows produced by dwarf 
galaxies, and since these dwarf galaxies are much more 
numerous than high-mass ones, they account for most of 
the IGM enrichment due to SNe. AGN-driven outflows 
are produced by massive galaxies, because these galax- 
ies are the ones that harbor AGNs. In this paper, we 
focus on SNe-driven outflows. AG N-driven outflows are 
presented in two separate papers (jGermain et al.l [20091 : 
iBarai et al.ll2010D . 

Numerical simulations and observations reveal that 
galactic outflows tend to be highly anisotropic. Sev- 
eral authors have performed SPH simulations of explo - 
sions inside isolated galaxies (|Mac Low fc Ferr ara 1999), 
or inside p rotogalaxies embeded in larger cosmological 
structures (jMartel fc Shapiro! l2001aHbh . These simula- 
tions reveal that outflows tend to be bipolar, with the 
energy and metal-enriched gas being channeled along the 
direction of least resistance. These results are supported 
by many observations that suggest that outflows are 
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anisotropic and even bipolar (e.g. iBland fc Tull 
Fab biano et al.lll990t|Shopbell fe Bland-Hawthor 
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also indirect support for anisotropic outflows from ob- 
servations of metal enrichment at high redshift: met- 
als are found in average-density regions far from known 
galaxies and a substantial scatter in metallicity is seen, 
which does not seem to b e purely a conseq uence of den- 
sity or galaxy proximity (|Pieri et all 120061) . Anisot ropic 
outflows travel larger distances than isotropic ones, and 
can more easily reach some low-density regions, while 
leaving others effectively pristine. 

Several numerical studies of the impact of galactic out- 
flows on the evolution of the IGM have been performed. 
The different approaches used fall into three categories. 
The hrst approach consists of describing the growth of 
large-scale structure and the formation of galaxies in the 
universe using either an N-body simulation or a semi- 
analytical method (e.g. Scannapicco & Broadhurst 2001, 
hereafter SBOl. lTheuns et al.ll200lUBertone et al.ll2005l ). 
and to combine it with an analyt ical model for describ - 
ing the evolution of the outflow ()Tegmark et al.l [l993f ). 
The model used assumes that outflows remain isotropic 
as they propagate through the IGM and the surrounding 
structures. 

The second approach uses a smoothed particle hydro- 
dynamics (SPH) algorithm to simulate both the forma- 
tion of large-scale structure and galaxies and the out- 
flows themselves. Outflows are generated by either de- 
positing additional thermal energy into SPH particles at 
the location of gala xies, to represent the energy gener- 
ated by SNe (e.g. iTheuns et aT1l2002bD . or by creat- 
ing a shell of SPH particles around galaxies, and giv- 
ing to th ese particles a large outward velocity compo- 
nent ( e.g. IScannapieco et alJteOOillSpringel fe Hernquistj 
120031 : lOppenheimer fe DaveH2006( ) . Outflows that start 
isotropic will become anisotropic as they propagate into 
a non-uniform external medium around galaxies. How- 
ever, unlike the analytical outflow model used in the first 
approach, SPH simulations have a limited resolution. If 
the anisotropy of the outflow is not caused by the den- 
sity distribution around the galaxies, but instead by the 
structure of the galaxies, the simulations will not be able 
to resolve it properly. 

The third approach consists of identifying galaxies in 
an output from an SPH simulation and calculating the 
propagation of outflows from these galaxies in N a differ- 
ent directions (|Aguirre etaD[2001). Since the resistance 
encountered by the outflows will be direction-dependent, 
outflows will start isotropic but then become anisotropic 
as the distance traveled by outflows will vary with direc- 
tion. In this approach, the matter located in the outflow 
slows down, but keeps expanding radially, while in the 
second approach described above, the matter located in 
outflows can be redirected in a different direction when 
it encounters resistance. But the most important limita- 
tion of this third approach is that it ignores the potential 
effect of feedback on the formation of galaxies, since the 
outflows are introduced a posteriori into an already com- 
pleted SPH simulation. 

To study anisotropic outflows in a cosmological con- 
text, including the effect of feedback, we have returned to 
the first approach. The analytical model used in these 
studies for describing the outflows do not suffer from lim- 



ited resolution, but a ssume isotropy. In a previous pa- 
per (jPieri et al.ll2007l hereafter Paper I), we presented a 
modified version of the outflow model, designed to de- 
scribe the evolution of anisotropic outflows. We then 
combined this model with the analytical Monte-Carlo 
method of SB01. We used a filtered Gaussian density 
field to determine the location, mass, and collapse red- 
shift of galaxies. Galaxies produce anisotropic outflows 
that propagate along the direction of least resistance. 
When an outflow hits a region destined to collapse by 
z = 2, it either strips it of its gas, preventing the forma- 
tion of a galaxy, or else enriches it in metals, potentially 
changing its cooling rate and the time necessary to turn it 
into a galaxy. We found that the anisotropy significantly 
enhanced the metal enrichment of average-to-low-density 
regions, as outflow tend to propagate into low-density re- 
gions, away from the cosmological structures (pancakes 
and filaments) where galaxies reside. 

The Monte Carlo approach used in Paper I provides a 
simple description of the formation and growth of struc- 
tures in the Universe, but it has a few drawbacks. First, 
the mass spectrum of halos is discrete, with only 10 differ- 
ent masses being allowed. Second, the peculiar velocity 
of halos is ignored; halos form at the comoving location 
of density peaks, and remain there throughout the course 
of the simulation. Hence, the clustering of halos and its 
potential consequences are poorly described. Third, the 
movement and accretion of gas on large-scales is ignored. 
Fourth, the treatment of halo destruction by mergers is 
quite simplistic. To address these four issues, we replace 
the semi-analytical approach of Paper I by a full N-body 
simulation of structure formation. 

This paper is set out as follows. In §2, we describe 
our numerical method and analytic outflow algorithm. 
Results are presented in §3, with and without photoion- 
ization suppression of galaxy formation. Implications are 
discussed in §4 and conclusions are presented in §5. 

2. THE NUMERICAL METHOD 

2.1. The Cosmological Simulations 

We perform a numerical simulation of dark mat- 
ter structure on cosmological scales by considering a 
ACDM model with present density parameter — 
0.268, baryon density parameter fi^n = 0.0441, cos- 
mological constant Ao = 0.732, Hubble constant Hq = 
70.4kms- 1 Mpc~ 1 (h = 0.704), primordial tilt n s = 
0.947, and CMB temperature Tomb = 2-725 cons istent 
with the results of WMAP3 4 (jSpergel et al.ll200l . We 
simulate structure formation inside a comoving cubic vol- 
ume of size Lbox = 15Mpc, with periodic boundary con- 
ditions. The simulation was perf ormed with a Particle- 
Parti cle /Particle-Mesh algorithm (jHocknev fe Eastwood! 
Il981f ) . using 384 3 equal-mass particles and a 1024 3 grid. 
The total mass in the computational volume is M to t = 
1.244 x 10 14 M Q , and the mass per particle is M part = 
2.197 x 1O 6 M . The length resolution is 5.86 kpc comov- 
ing. 

2.2. Halos, Merger Trees, and Galaxy Formation 

The simulation produces dumps containing the posi- 
tions and velocities of the particles. The dumps are 

4 http: / /lambda. gsfc.nasa.gov/product/map/dr2/params/lcdm_all.cfm 
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regularly spaced in time by a time interval At = 2.5 x 
10 7 years, which gives 131 dumps between redshifts z = 
24 and z = 2. In each dump, we identify halos using the 
friends-of- friends algorithm (|Davis et al.l Il985t ). with a 
linking length equal to 0.2 times the mean particle spac- 
ing. As in Paper I, we set the minimum mass of halos 
at M m j n = 7.61 x 10 7 M Q , c orresponding to 34 pa rticles. 
As we found in Paper I and lPieri fe Martell (|2007t) . halos 
below this mass do not form galaxies with mature stellar 
populations either because of long cooling times or pho- 
toionization heating. We do not consider the impact of 
mini-halos or Population III stars here, we assume that 
it is negligible aside from pre-enriching halos to a level 
where mature stellar populations can arise. 

By combining the catalogs of halos, we can find the 
ancestry of each halo. Consider a particular halo H™ 
found in dump n. We identify the particles contained in 
this halo, and then examine the previous dump, n — 1, 
to find where these particles were. If 50% or more of the 
particles located in a halo in dump n — 1 end up in 

halo H™ in dump n, then halo -ff™ -1 is considered to be a 
progenitor of halo H™. Then, we have three possibilities, 
depending on the number of progenitors: 

• Case A: Halo Hf has no progenitor. This halo 
was formed by the monolithic collapse of a region 
containing field particles, clumps (halos containing 
fewer than 34 particles), or pieces extracted from 
other halos that amounted to less than 50% of those 
halos. In this case, we consider that a new halo is 
born, the birth redshift being the redshift corre- 
sponding to dump n. Halo H™ is a "leaf" in a 
merger tree. 

• Case B: Halo H™ has one progenitor. In this case, 
the halo already existed in dump n — 1. Between 
the two dumps, it might have grown by accretion of 
field particles or small clumps, or might have lost 
mass by evaporation. 

• Case C: Halo H™ has several progenitors. In this 
case, the merger of the progenitors have formed 
a new halo, and the progenitors no longer exist. 
There is an exception: if one progenitor, say halo 
flJT 1 , provides 90% or more of the mass of halo 
H™, we do not consider this to be a merger. Halos 

and iJ™ are actually the same halo, which 
grew by accreting (not merging with) smaller halos. 
In this case, halo if™ is not treated as a new halo, 
but as a preexisting one. 

Figure Q] illustrates the various cases. Halo b, which 
had no progenitor, formed by monolithic collapse (also, 
halo d, and halo a in dump n — 3). Halo a grows by 
accretion between dumps n — 3 and n — 2. Halo c formed 
by the merger of halos a and b. Halo d is accreted onto 
halo c. This is not regarded as a merger because 90% 
or more of the mass of halo c in dump n came from one 
single progenitor. 

Once a halo has formed, if its gas is not being kept 
hot by the UV background, it will lose energy by radia- 
tive cooling. We implement this photoionization heat- 
ing in manner descr ibed by the fiducial approach of 
iPieri fe Martell (|2007l) . This is manifest in the simula- 
tion as a redshift-dependent minimum halo mass for gas 
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Fig. 1. — Merger tree illustrating the various cases: monolithic 
collapse (Case A), growth (Case B), merger (Case C), and accretion 
(Case C, exception). The areas of the circles are proportional to the 
masses of the halos. Lowercase labels identify halos, with identical 
labels indicating halos that existed in the previous dump. 



cooling using the mass for which half the halo ga s is ex- 
pected to cool (taken from iDiikstra et al.l 120041 ). This 
corresponds to a halo circular velocity of 

"1/2= [89.4 ±2.4- (6.3 ±0.3)z con ]kms" 1 , (1) 

where z co \\ is the collapse redshift of the halo. We conser- 
vatively apply this criterion by requiring that halos must 
have reached turnaround after z — 6 for suppression of 
cooling to occur. We investigate the impact of photoion- 
ization heating in iJ3] by showing results, first without, 
then with this criterion. 

Where photoionization suppression does not occur, 
halo gas cools resulting in galaxy formation and outflow 
production. As in SB01 and Paper I, we calculate the 
virial temperature of each hal o from its mass a nd collapse 
redshift using the method of Eke et al. (1996), and then 
we calculate the coolin g time using the cooling model of 
IWhite fc Frenkl (|1991f ). We assume that the halo turns 
into a galaxy after the cooling time has elapsed. As it 
turns out, for the range of halo masses and collapse red- 
shifts considered in our simulations, the cooling times 
are almost always very long or very short, that is, ei- 
ther longer than the age of the universe or shorter than 
the time interval At = 2.5 x 10 7 years between dumps. 5 
Hence, most collapsed halos either become galaxies al- 
most immediately, or never become galaxies. We assume 
that once a galaxy has formed, the outflow starts imme- 
diately. This is a valid assumption, since the lifespan of 
the progenitors of Type II SNe are short compared to the 
time resolution of our simulations. 

Mergers and accretion of halos onto more massive ha- 
los bring an end to the fueling of outflows from hosted 
galaxies, since these galaxies no longer exist as indepen- 

5 This is caused by the rapid drop in cooling rate at temperatures 
T < 10 4 K. 
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dent systems. The effect is to revert any outflows into 
a post-supernova phase (see the following section) and 
begin the process of describing the properties of the new 
merged galaxy as a source of outflow fueling. 

2.3. The Outflows 

2.3.1. Direction of Least Resistance 

To determine the direction of least resistance out of 
a halo, we first determine the orientations and lengths 
of the semimajor axe s of the hal o , usin g the standard 
method described by iDavis etldl (|1985f ) . We calculate 
the quadrupole tensor Q of the halo, 6 



Qij = y^mfc(r fc )i(r fc )j 



h3 =x,y,z: 



(2) 



where and are the mass of particle k and its posi- 
tion relative to the center of mass, respectively. We then 
diagonalize this tensor. The semimajor axes oi, a,2, a% 
of the halo are related to the eigenvalues Qi, Q2, Q3 by 
Q>% = (SQi/M) 1 / 2 , where M is the mass of the halo, and 
the eigenvectors give us the direction of these axes. We 
then take the direction of least resistance as being along 
the shortest axis. 

Wc have to make sure that the number of particles per 
halo is sufficient for this method to work. In Paper I, 
the smallest halos we considered had a mass of 7.61 x 
10 7 M Q . Here a halo of such mass would be made of 34 
particles, which is sufficient to determine the direction of 
the shortest axis. 

Figure [5] shows a region of size 2.1 Mpc and thickness 
0.25 Mpc at redshift z = 4 in the computational volume. 
We focus on a region containing a large-scale cosmolog- 
ical filament. The red dots represent the P 3 M particles, 
the black dots indicate the positions of clusters identi- 
fied by the friend-of-friend algorithm, and the bars in- 
dicate the direction of least resistance. These bars have 
equal length in 3D, but since they are seen in projection, 
they appear to have different lengths. There is clearly 
an overall tendency of the directions of least resistance 
to be aligned with each other, and to be perpendicu- 
lar to the filament. Outflows propagating along these 
directions will transport energy and metal-enriched gas 
preferentially in low-density regions. 

2.3.2. Expansion 

The dynamical equations describing the evolution of 
the outflows are presented in full detail i n Paper I. They 
are ba sed on the original treatment of [Tcgmark et al.l 
(|1993t ). SBOl, and IScannapieco et al.l (|2002j), and have 
been generalized to the case of anisotropic outflows. Here 
we present a review of the method. 

The injection of thermal energy by SNe produces an 
outflow of radius R, which consists of a dense shell of 
thickness RS containing a cavity. A fraction 1 — f mi of 
the gas is piled up in the shell, while a fraction, f m , of the 
gas is distributed inside the cavity. We normally assume 
that most of the gas is located inside a thin shell (S <C 1, 
fm <C 1)- This is called the thin-shell approximation. 

The evolution of an outflow of radius R and opening 
angle a expanding out of a halo of mass M is described 

6 Incorrectly called "inertia tensor" by many authors. 




Fig. 2. — Region of size 2.1 Mpc and thickness 0.25 Mpc centered 
on a large-scale filament, at redshift 2 = 4. Red dots: P 3 M par- 
ticles; Black dots: halos. The lines indicate the directions of least 
resistance. 

by the following system of equations, 
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where a dot represents a time derivative, f2, fif,, and H 
are the total density parameter, baryon density param- 
eter, and Hubble parameter at time, t, respectively. L 
is the luminosity (see below), p is the thermal pressure 
resulting from this luminosity, and p G xt is the external 
pressure of the IGM. a is the opening angle of the out- 
flow, which is 180° for isotropic outflows. 

The external pressure p ext depends upon the density 
and temperature of the IGM. As in Paper I, we assume a 
photoheated IGM made of ionized hydrogen and singly- 
ionized helium, wit h a fixed temperature Tiqm = 10 4 K 
(|Madau et al.| [2001). and an IGM density p IGM equal to 
the mean baryon density, pb- The external pressure is 
then given by 

PbkTiQM 3£l fc ,o#ofcTiGM(l + z) 3 

Pextiz) = = — ■ — Vt; ' ( 5 ) 

where z is the redshift and /1 = 0.611 is the mean molec- 
ular mass. 

The luminosity L is the rate of energy deposition or 
dissipation within the outflow, and is given by 

L(t) = isN — icomp i (6) 

where LgN is the total luminosity of the supernovae re- 
sponsible for generating the outflow, and L comp repre- 
sents cooling due to Compton drag against CMB pho- 
tons (other contributions are neglected, as in Paper I). 
The supernovae luminosity, for a galaxy forming in a halo 
of mass M, is given by 



L 



SN 
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where /» is the star formation efficiency, A/* = 
f*MQb t o/Qo is the total mass in stars formed during the 
starburst, f w is the fraction of the total energy released 
that goes into the outflow, M req is the mass of stars re- 
quired to form one SN, and Eq = 10 51 ergs is the energy 
released by each of these SNe. As in Paper I, we take the 
values /* = 0.1, A/ rcq = 89.7M^ (deriv ed using a bro- 
ken power-law IMF from [Kroupa 2 00lf). and the mass- 
depen dent expression for / w given by IScannapieco et al.l 

(EqqI). 

The Compton luminosity for a mixture of ionized hy- 
drogen and singly-ionized helium is given by 



L, 



comp 
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where <r t is the Thompson cross section, and T 7 o is the 
present CMB temperature. We used T 1 = T 70 (l + z), 
which is valid over the range of redshifts we consider. 

The expansion of the outflow is initially driven by the 
supernovae luminosity. After a time tburst — 5 x 10 7 yr, 
the supernovae turn off, and the outflow enters the "post- 
supernova phase." The pressure inside the outflow keeps 
driving the expansion, but this pressure drops since there 
is no energy input from supernovae. Eventually, the pres- 
sure will drop down to the level of the external IGM pres- 
sure. At that point, we assume that the expansion of the 
outflow will simply follow Hubble expansion. 

2.3.3. Ram Pressure Stripping 

As outflows propagate into the IGM, the expanding 
shells of swept-up gas may eventually hit halos. Fol- 
lowing SB01 and Paper I, we assume that if an already 
formed galaxy is hit by an outflow, the cross section is 
too small for the impact to have any significant effect. If 
a halo that has not yet collapsed and formed a galaxy is 
hit, two things might happen. Either the ram pressure 
of the outflow will strip the halo of its baryonic content, 
preventing the formation of a galaxy, or else the halo will 
be enriched in metals by the outflow. The condition for 
stripping is 



I 2 \ 

1 M v > M b v c 



4R 2 J 



(9) 



where M a and R the mass and radius of the shell, re- 
spectively, v is the outflow velocity, I and M& are the 
radius and baryonic mass of the halo being hit, and v esc 
is the escape velocity from the halo. In Paper I, we used 
the spherical collapse model to estimate the radius I at 
the time of the hit. This model predicts that halos ex- 
pand, eventually turn around, and collapse to a point. 
This was straight forward in the Monte Carlo approach 
of Paper I. However, in the numerical approach used in 
this paper, we identify collapsed halos using the friend- 
of-friend algorithm, and these halos obviously have a fi- 
nite size. Thus, the basic spherical model needs to be 
modified to account for the fact that real halos do not 
collapse down to a point, but instead contract by a factor 
of order 2 after turnaround, and reach virial equilibrium. 
We describe our modified spherical collapse model in Ap- 
pendix B. 

When the criterion for stripping is met, the halo is 
labeled as being stripped, and it will not form a galaxy. 



If the halo is not stripped, the outflow will deposit metals 
into the halo. 

2.3.4. Metal Enrichment of Halos 

We assume that halos start with metallicity of 
[Fe/H] = —3, which is negligible for the purposes of 
calculating the halo cooling time. Once galaxies are 
formed, metals are pro duced at rate of 2M Q per SN 
(|Nagataki fc Sato| [l998'). Hence the mass of metals in 
the outflow is 



M z = f & 



(10) 



where / esc is the fraction of ISM gas blown out with the 
outflow. We use the va lue / eS c = -5 tak en from the 
numerical simulations of iMori et al.1 (|2002[ ). This mass 
of metals is distributed evenly throughout the volume of 
the outflow. 

When an outflow strikes a halo and does not strip it, 
it modifies the metal content of the halo by depositing a 
fraction of its metals, /de P Overlap /Kutflow, where / dep is 
a mass deposition efficiency, which we set at /d cp = 0.9, 
and Overlap is the volume of the outflow that overlaps 
with the halo. For details, we refer the reader to Paper I. 

The radiative cooling rate increases with metallicity. 
So when a halo is enriched by metals, we recompute its 
cooling time, which determines when this halo will form a 
galaxy. This is the potential to "bring to life" dead halos. 
A halo that cannot become a galaxy because its cooling 
time exceeds the age of the universe might form a galaxy 
after all, if metal enrichment reduces the cooling time. It 
turns out to be a rare occurrence. In all the simulations 
we performed only several galaxies were formed this way, 
out of 20 000 - 60 000 galaxies. 

2.4. Metal Enrichment of the IGM 

Stellar evolution in galaxies produces metals, that are 
then carried into the IGM by outflows. Of particular in- 
terest is the volume filling factor, that is, the fraction of 
the IGM, by volume, that has been enriched in metals 
by outflows. To calculate the volume filling factor, we 
cannot simply add up the final volumes of the outflows, 
and divide by the volume of the computational box. In- 
tergalactic gas enriched by outflows will move with time 
as structures grow, and therefore regions that were never 
hit by outflows might end up containing metals. 

We take this effect into account, by using a dy- 
na mic particle en r ichme nt scheme that was developed 
bv lGermain et al.l f|2009h ■ By combining the dumps pro- 
duced by the P 3 M code, which contains the position of 
the particles, with the positions, orientations, and radii 
of the outflows, we can find which particles are being 
hit by outflows, and at what redshift. We can then use a 
Smoothed Particle Hydrodynamics technique to estimate 
the volume filling factor at any redshift. A smoothing 
length h is ascribed to each particle. We calculate ft. it- 
eratively by requiring that each particle has between 60 
and 100 neighbors within a distance 1.7ft.. We then treat 
each particle as an extended sphere of radius 1.7ft. over 
which it is considered to be spread. We then divide the 
computational volume into N = 256 3 cubic cells. The 
cells that are covered by one or more enriched particles 
are then considered enriched. Notice that simply count- 
ing the number of enriched particles in each cell would 
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Fig. 3. — Volume filling factor versus redshift, with various open- 
ing angles, for simulations without photoionization suppression of 
galaxy formation. 

fail in low-density regions, where the local particle spac- 
ing exceeds the cell size 

Additionally, giving a smoothing length to each P 3 M 
particle enables us to easily calculate the gas density in 
the center of each cell, using the standard SPH equation, 

<°( r g) = ^ m iW(r g - Tj, hj) , (11) 

i 

where r,-, m„ and hi are the position, mass, and smooth- 
ing length of particle i, r g is the pos ition of a grid po int, 
and W is the smoothing kernel (see I Monaghanl 1 1 9921 for 
a detailed description of SPH). We will use this in § 3.2 
below. 

3. RESULTS 

We performed a series of ten simulations: five with 
our galaxy formation scheme excluding suppression by 
photoionization heating, to produce results that can be 
compared with our results in Pa per I, and a further fiv e 
with this suppression as used in iPieri fc Martell (|2007f ). 
In each case we simulate outflows with opening angles 
a = 60°, 90°, 120°, 150°, and 180°. This latter case 
corresponds to isotropic outflows. 

3.1. Without Photoionization Suppresion 
3.1.1. Volume Filling Factor 

Figure [3] shows the evolution of the volume filling fac- 
tor for the five opening angles considered. Enrichment 
starts at redshift z ~ 12 when the first galaxies form, and 
steadily increase with time, to reach 8 — 28% by redshift 
z = 2. The fraction increases with decreasing angle, at 
all redshifts. This differs from the results of Paper I, and 
clearly shows the importance of an accurate description 
of clustering. In Paper I, we found that the fraction was 
nearly constant, near 14%, for angles between a = 180° 
and a = 100°, and then dropped down to 11% at a = 60° 
(Fig. 8, bottom left panel, in Paper I). 

Three distinct effects determined the overall volume 
filling factor in Paper I. These effects are still present 
here but their relative impact is changed. First, out- 
flows are much more likely to overlap when the sources 
are highly clustered and this effect is diminished where 
outflows are anisotropic and aligned. Second, clustering 
makes it more likely that an expanding outflow will hit 
a pre-galactic collapsing halo, possibly resulting in gas 



TABLE 1 

Statistics of Enrichment and Hits Without 
Photoionization Suppresion 



a 


^outflows 


V.F.F. 


N ut 


-^stripped 


-^outflows 


60° 


0.388 


0.277 


159 957 


145 188 


63 015 


90° 


0.381 


0.206 


192115 


168 365 


52 354 


120° 


0.369 


0.155 


213 846 


178 304 


44 399 


150° 


0.367 


0.109 


231 058 


178 781 


38 707 


180° 


0.369 


0.081 


234 365 


165 353 


34 324 



stripping and preventing a galaxy from forming and pro- 
ducing its own outflow. Again the impact of this effect 
is diminished where outflows are anisotropic. Third, the 
volume of each outflow decreases with decreasing angle, 
the increase in radius not being enough to compensate 
for the smaller angle. The first two effects tend to in- 
crease the volume filling factor for anisotropic outflows 
and are sensitive to the level of galaxy clustering. The 
third effect tends to decrease the overall volume filling 
factor and is insensitive to the degree of clustering. Our 
results here show a more accurate and stronger degree 
of halo clustering, thus explaining the more marked im- 
pact of the first two effects with respect to the third. 
This also explains the small fractions we obtain in the 
isotropic case here (8% at z = 2) by comparison to the 
equivalent result in Paper I. 

In Table [TJ we list, for each opening angle, the total 
volume of the outflows, Voutfiows (sum of all the volumes, 
ignoring overlap, in units of the box volume), the vol- 
ume filling factor (V.F.F.) the number of halos hit and 
stripped (A^itj -Stripped, respectively), and the number 
of outflows, ^outflows, all at z = 2. Comparing Outflow 
with V.F.F. in Table [TJ we clearly see that the effect of 
overlap is mild for small opening angles, but becomes 
very important for large ones. As the opening angle in- 
creases, there are more halos hit, more halos stripped, 
and as a result fewer galaxies are formed, producing 
fewer outflows. Going from a — 60° to a = 180°, the 
number of outflows drops by a factor of 2. Despite this, 
though, their total volumes are comparable (0.388 vs. 
0.369) since the volume per outflow is greater for out- 
flows with larger opening angles. 

In addition to the three effects described above we are 
also able to consider a fourth effect given the dynami- 
cal nature of our new simualtions. The matter enriched 
with metals moves as the simulation proceeds, and this 
motion modifies the volume filling factor. An isotropic 
outflow produced by a galaxy enriches matter that is lo- 
cated near the galaxy, and this matter is likely to later 
contract as it accretes over the cosmological structure 
hosting the galaxy. Highly anisotropic outflows travel 
larger distances, and might enrich matter located far 
from cosmological structures, and potentially expanding 
away from them. 

Notice that the number of galaxies formed is signifi- 
cantly larger than in Paper I (~ 20 000 in a cosmological 
volume that was 40% larger). As we showed in Paper I, 
the Monte Carlo method underestimates the number of 
galaxies compared to either an N-body simulation or the 
Press-Schechter approximation (see Fig. 4 in Paper I). 
This provides additional motivation for performing a full 
cosmological simulation. 
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3.1.2. Distribution of Metal-Enriched Gas 

Figures HHHl show slices of the computational volume 
at redshifts z = 4.98, 4.00, 3.00, and 2.00, for opening 
angles a = 180°, 120°, and 60°, respectively. At high 
rcdshift, metals are located in high-density regions where 
the sources are located. As time goes on, they eventually 
reach lower-density region, but never travel far enough 
to reach the deepest voids. Instead, even at z = 2, all 
metals are found near dense structures like clusters and 
filaments. 

The main effect of narrowing the opening angle is to 
increase the volume over which metals are dispersed, 
as Figure [3] showed. Comparing Figures HHH it ap- 
pears that both high- and moderately-low-density re- 
gions are affected. It is not obvious from these fig- 
ures that anisotropic outflows enrich predominantly low- 
density regions. To investigate this issue, we calculated 
the density at each of the N cells on the grid, using equa- 
tion (TiTl) . We then binned these cells according to den- 
sity, and counted in each bin the number of cells, N p , and 
the number of enriched cells, N' p . From this, we calcu- 
lated statistics of metal enrichment vs. density. The re- 
sults are shown in Figure [7J The top panel shows N' /N p ; 
this is effectively the probability of enriching a systems 
of a given density. The galaxies producing the outflows 
are located in high-density regions, hence these regions 
are favored over low-density ones, even with anisotropic 
outflows (such outflows will still enrich the gas located 
near the galaxies). The enrichment reaches 80% or more 
in high-density regions. The extent of enrichment in- 
creases with decreasing opening angle, at all densities, 
although this effect is most significant at lower densities. 
Going from a = 180° to a = 60°, N' p /N p increases by 
a factor of 3.90 at log[p/p] = — 1, 3.53 at log[p/p] = 0, 
1.75 at log[p/p] = 1, 1.24 at log[p/p] = 2, and 1.13 at 
log[p/p] =3. 

The bottom panel shows the number of grid points en- 
riched at a given overdensity as a fraction of the total 
number of cells, N . There are very few cells at densi- 
ties below log [p/p] = —1.2, but about ~ 10% of these 
cells get enriched, explaining why N'/N becomes negli- 
gible while N' p /N p does not. The curves are very similar 
for all opening angles, except for the overall amplitude. 
Most enriched cells are located at low-densities in the 
range log [p/p] = — 1 to 0, even though enrichment is less 
efficient in these regions compared to high-density ones. 
Even in the case of isotropic outflows most of metals are 
in underdense regions when considered by volume. 

As the opening angle is reduced, outflows travel pre- 
dominantly into low-density regions, and the number of 
outflows increases, because less stripping occurs. These 
two effects are acting in the same direction in low-density 
regions, while they are competing in high-density re- 
gions. To separate the two effects, we replotted in Fig- 
ure [8] the same results as in Figure but rescaled by 
the number of outflows. That is, we have multiplied 
the number of enriched cells N' by the factor /outflow = 

A r outflow(180°)/A^outflow S (a), where AUtftow can be read 
from the last column of Table [T] The top panel of Fig- 
ure |8] is strikingly similar to the top panel of Figure 9 
in Paper I. As the opening angle decreases, low-density 
regions are being enriched more, at the expense of high- 



TABLE 2 

Statistics of Enrichment and Hits With 
Photoionization Suppresion 



a 


^outflows 


V.F.F. 


N ut 


-^stripped 


-^outflows 


60° 


0.301 


0.189 


84 599 


80 082 


49 276 


90° 


0.280 


0.119 


107 242 


100 410 


39 994 


120° 


0.255 


0.079 


120 593 


110 588 


32 773 


150° 


0.238 


0.065 


126114 


112 347 


27 339 


180° 


0.223 


0.047 


120 717 


103 479 


22 992 



density regions. However, the crossover happens at a 
higher density, p/p « 10, compared to p/p ~ 2 if we 
make a similar correction to /outflow for Paper I. We 
attribute this difference to our dynamical treatment of 
metal enrichment. Enriched gas located in low-density 
regions can eventually move into higher-density regions 
as the system evolves. 

3.2. With Photoionization Suppression 

Figure [9] shows the evolution of the volume filling fac- 
tor with the inclusion of photoionization suppression of 
galaxy formation as described in ^2.21 By comparison 
with Figure [3j one can see that the volume filling factor 
remains unaffected until z w 4 and after this the extent 
of mechanical feedback is diminished. The volume fill- 
ing factor continues to increase until z w 3 after which 
it actually begins to fall until the end of our simulations 
at z = 2. This is a result of the dynamics included in 
our simulation, whereby gravitational attraction carries 
metal-enriched regions towards dense structures faster 
than outflows can carry them away from these structures. 
The final statistics of the suite of runs with photoioniza- 
tion suppression are shown in Table [2] The same nota- 
tion as Table [T] is used and as before increasing opening 
angles leads to more outflows hitting halos, more strip- 
ping and so fewer galaxies (and so outflows). 

As one can see by comparing Table [2] with Table [TJ 
the introduction of photoionization suppression has the 
weakest impact for our smallest opening angle: the vol- 
ume filling factor is 32% lower due photoionization heat- 
ing in the for a — 60°. However, the most pronounced 
impact occurs one of our intermediate models, a = 120° 
where the volume filling factor is 49% lower. The number 
of halos hit is also non-linear with opening angle, peaking 
at a = 150°. It is also notable that the summed volume 
of all outflows Voutflows shows a clear decline with in- 
creasing opening angle, indicating that the fall in galaxy 
numbers with increasing opening angle dominates over 
the impact of the change in individual outflow volumes. 
As in the case without the inclusion photoionization sup- 
pression, the true model volume fraction shows a stronger 
opening angle dependence and this is a sign of reduced 
outflow overlap for more anisotropic outflows. 

The density dependence of the volume filling factor is 
shown in Figure [T0l These results are broadly similar 
to those for the suite of simulations without photoion- 
ization suppression (Figure [7]). The values of N' p /N are 
reduced in underdense regions (log [p/p] < 0), but almost 
unaffected in very high-density regions (log [p/p] > 1). In 
these high-density regions, outflows strongly overlap, and 
the removal of some of them by photoionization suppres- 
sion does not prevent cells from being enriched. 



Fig. 4. — Slice (of comoving size 15Mpc X 15Mpc, and comoving thickness O.lMpc) of the computational volume for an opening angle 
a = 180°, showing the evolution of the distribution of metals, for simulations without photoionization suppression of galaxy formation. 
The areas shown in green are enriched. The black dots represent the P 3 M particles, and show the large-scale structures. 



Fig. 5. — Same as Figure[4]for an opening angle a = 120°. 



Fig. 6. — Same as Figures [4] and \E\ for an opening angle a = 60°. 
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Fig. 7. — The number of enriched cells N' p in the simulation 
volume at z = 2 as a function of IGM density with various open- 
ing angles, for simulations without photoionization suppression of 
galaxy formation. Top panel: N' p as a fraction of the number of 
cells at this density, N p . Bottom panel: N' p as a function of the 
total number of cells, N . 
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Fig. 8. — Same as Figure [7] except that the number of 
enriched cells has been multiplied by the factor /outflow ( a ) = 
A r outflow(180 o )/Af O utflows(a), where JV outfiows is given in Table [l] 
This isolates the impact of anisotropy on the regions enriched by 
individual outflows. 

4. DISCUSSION 

In Paper I, we described the evolution of anisotropic 
outflows and the enrichment of the IGM using a sim- 
ple Monte Carlo method that was originally introduced 
by SB01. In this paper, we used cosmological N-body 
simulations to simulate the formation and evolution of 
large-scale structures in a more realistic way. The re- 
sults obtained with both methods show interesting dif- 
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Fig. 9. — Volume filling factor versus redshift, with various 
opening angles, for simulations with photoionization suppression 
of galaxy formation. 
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Fig. 10. — The number of enriched cells N' p in the simulation 
volume at z = 2 as a function of IGM density with various opening 
angles, for simulations with photoionization suppression of galaxy 
formation. Top panel: N' p as a fraction of the number of cells 
at this density, N p . Bottom panel: N' p as a function of the total 
number of cells, TV. 

ferences. In the Monte Carlo approach, the volume filling 
factors were largely independent of opening angle (Fig. 8, 
bottom-left panel of Paper I) . Here were find a clear de- 
pendence on outflow anisotropy: the volume filling fac- 
tor increases with decreasing opening angle, as Figures[3] 
and [9] showed. 

In Paper I, we generated a Gaussian density field at 
high redshift, filtered it at various mass scales, and used 
the spherical collapse model to predict the location and 
collapse redshift of halos. Hence, the location of halos 
was predetermined by the initial conditions, and these 
locations did not change with time. Also, the density 
p was determined by a nonlinear mapping of the initial 
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density field. The density p(r) at a comoving position r 
was thus entirely determined by the initial density at that 
location. This approach has the advantage of simplicity, 
but does not constitute a full treatment of clustering. 

This can have important consequences. Galaxies tend 
to be located inside larger cosmological structures, like 
filaments or pancakes. If several galaxies are located in 
a common structure, they can affect each other in two 
ways. First, the outflows might overlap, and second, if 
an outflow hits a halo that has not yet formed a galaxy, 
ram pressure stripping might prevent the galaxy from 
forming. These two effects are reduced when the outflows 
are anisotropic, because such outflows travel away from 
the cosmological structures, into low-density regions, and 
tend to avoid other galaxies and other outflows. 

These effects were present in the Monte Carlo simu- 
lations (Fig. 8, bottom right panel of Paper I), how- 
ever, they were not dominant. Instead they either ap- 
proximately balanced the decline in individual outflow 
volumes for increasingly anisotropic outflows (down to 
opening angles of a — 100°) or were dominated by them 
(opening angles from a = 100° to a = 60°). The strength 
of clustering is much more apparent in the cosmologi- 
cal N-body simulations, and these effects now dominate. 
The higher level of clustering in the simulations here give 
greater weight to the impact of outflows that travel pref- 
erentially into low-density regions. Most results which 
are sensitive to this effect are all stronger here: the re- 
duced outflow overlap, the reduced stripping of halos, 
the increased extent of enrichment in low-density regions. 
There is one aspect in which these results show a weaker 
signal of anisotropy than Paper I: the enrichment of over- 
dense systems. In Paper I the fact that outflows travel 
into low-density regions resulted in a reduction in the 
physical extent of outflows in overdense regions, here 
however, anisotropic outflows increase the extent of out- 
flows for densities all the way up to p/p « 10. This a 
consequence of the dynamic nature of our N-body simu- 
lations. At the time of launching the outflow may expand 
into underdense regions but as the simulation progresses 
the gas enriched by the outflow is re-accreted onto over- 
dense structures. 

The impact of this accretion is also notable in the con- 
text of results from lOppenheimer et ail (|2009[ ) . They 
find that understanding the accretion of outflow gas back 
onto galaxies plays an important role in galaxy assembly 
and the stellar mass function of galaxies by z = 0. We 
find accretion of outflow gas onto large-scale structures 
plays an important role at z > 2 and it seems likely that 
much of this gas continues its infall onto galaxies and 
fuels star formation. 

Consider a galaxy producing an outflow that carries 
metals into low-density regions. As Figures 0H6] show, 
outflows never travel very far from their sources. Thus, 
the gas enriched by outflows is likely to be gravitationally 
bound to the structure hosting the galaxy. As that gas 
accretes onto that structure, it carries metals from a low- 
to a high-densi t y reg ion. This effect has been found by 
I Germain etal] ([2009) in a different context (anisotropic 
outflows powered by AGNs). As these authors argue, 
even though anisotropic outflows deposit metals predom- 
inantly in low-density regions, the subsequent evolution 
of the large-scale structure will tend to partly wash out 
that effect. The only exception might occur when out- 



flows travel very large distances, reaching gas inside deep 
voids that is not destined to accrete onto any structure 
within a Hubble time. This hap pens with AGN-powered 
outflows (|Germain et al.l [20091 ). but SNe-powered out- 
flows do not have enough energy to reach such distances. 

The simulations presented in this paper use the same 
outflow model as the ones presented in Paper I. Only 
the method used for generating the population of halos 
(Monte Carlo vs. N-body simulations) and for calculat- 
ing the enrichment of the IGM (static vs. dynamical) 
are different. The outflow model contains several pa- 
rameters, such as /*, / w , /esc, and M roq . We use the 
same values as in Paper I, and we refer the reader to 
that paper for a discussion on the limits on these various 
parameters. 

5. SUMMARY AND CONCLUSION 

We have combined cosmological N-body simulations 
with a friend-of-friend algorithm for identifying halos 
and merger events, and an anisotropic outflow model, to 
study the metal-enrichment of the IGM by SNe-powered 
galactic outflows, in a ACDM universe. Outflows are 
modeled as bipolar cones traveling along the direction 
of least resistance, which we take as the direction along 
which the density drops the fastest. We performed five 
simulations with outflow opening angles ranging from 60° 
to 180°, the latter case corresponding to isotropic out- 
flows. Each simulation was stopped at redshift z — 2. 
For each simulation, we tracked the deposition of met- 
als in the IGM by outflows, and the subsequent motion 
of the metal-enriched gas. We also included the effect 
of ram pressure stripping and metal enrichment of pre- 
galactic collapsing halos that are hit by outflows. We 
did this with and without the inclusion of photoioniza- 
tion suppression of galaxy formation resulting in a total 
of ten simulations. Our main results are the following. 

• Anisotropic outflows travel predominantly into 
low-density regions, away from cosmological struc- 
tures (filament and pancakes) that host most galax- 
ies. As a result, anisotropic outflows are less likely 
to overlap than isotropic ones, and are less likely 
to inhibit galaxy formation by stripping collapsing 
halos of their gas. The combination of these effect 
results in an increase of the volume filling factor of 
the IGM with decreasing opening angle, from 8% 
with isotropic outflows, up to 28% with opening 
angle of 60°. 

• The majority of collapsing halos hit by outflows 
end up being stripped of their gas, preventing 
the formation of a galaxies. The total number of 
galaxies formed drops by a factor of 2 when go- 
ing from opening angle a — 60° to 180°. Halos 
not stripped are enriched in metals, with negligible 
consequences for the formation time of the galaxy. 
In particular, fewer than one in 10 3 - 10 4 galax- 
ies that formed would not have formed without the 
hit, because their cooling time exceeded the Hubble 
time. 



• The volume filling factor remain small (28% or 
less) because outflows do not travel large distances. 
They remain fairly close to the structures where 
they originate, never reaching the center of the 
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deep cosmological voids. The regions most effi- 
ciently enriched (up to 90% by volume) are the 
high-density regions, where the galaxies producing 
the outflows are located. The enrichment of low- 
density regions is not as efficient. However, by vol- 
ume, most metals are located in these low-density 
region, because their combined volume greatly ex- 
ceeds the combined volume of high-density regions. 
Indeed, most of the enriched volume is in regions 
whose density is between O.lp and p. 

• As outflows become more anisotropic (smaller 
opening angle a), the extent of enrichment in- 
creases at all densities, low and high. The effect 
is most important at low densities, where the effi- 
ciency increases by factors of several. At very-high 
densities, p — lOOOp, regions are already 80% en- 
riched by volume by isotropic outflows, so there is 
little room for an additional increase. 

• Overall, anisotropic outflows do not strongly fa- 
vor low-density regions at the expense high-density 
ones, unlike what was found in Paper I. Anisotropic 
outflows are more numerous, resulting in higher 
probability of enrichment at all densities. When 
this effect is factored-out, we find that the prob- 
ability that a galaxy will enrich systems at densi- 
ties up to 10/O is higher for increasingly anisotropic 
outflows. This results from evolution of the metal- 
enriched IGM. Enriched gas tends to be located rel- 



atively near the large-scale structures hosting the 
galaxies responsible for this enrichment. As that 
enriched gas accretes onto these structures, the ef- 
fect of anisotropy is partly washed out. 

• Photoionization suppression of low-mass galaxy 
formation leads to a number of modifications to 
our enrichment predictions. Since many of these 
low mass galaxies would otherwise have formed in 
or near voids of our density distribution, their sup- 
pression leads to a reduction in the enrichment of 
low-density regions. Also a surprising, but not im- 
plausible, consequence of this prescription is a fall 
in volume fraction of enriched regions after z = 3. 
This is a consequence of suppression of much late 
galaxy formation combined with the accretion of 
enriched structures back onto high-density, low- 
volume regions. 
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APPENDIX 

MERGER TREE: DEALING WITH SPLITTERS 

While building merger trees, we discovered that some clusters of P 3 M particles were actually splitting, that is, 
breaking up into two or more components. Such splitting might be real. A small cluster located in the vicinity of a 
larger one might get stretched by the tidal field to the point where it breaks up, or it might get so elongated that 
the FOF algorithm identifies it as several clusters. This can also be a purely numerical effect, which could occur in 
two different situations: (1) Two clusters that pass near each others without merging come briefly into contact. If 
the FOF algorithm is applied while the clusters are in contact, they will be identified as one single cluster, that then 
splits. (2) When two clusters approach each others on a collision course, they will be identified as one single cluster 
as soon as one particle from the first cluster and one particle from the other cluster are within a linking length of 
one another. Then, because of the internal motion of the particles inside clusters, the two clusters might lose and 
re-establish contact several times during the following timesteps, before they finally merge for good. These effects are 
artifacts of the FOF algorithm, and usually go unnoticed unless one uses a very fine time resolution when building 
merger trees. In this paper, we build cluster catalogs spaced by time intervals of 2.5 x 10 7 years, which gives us merger 
trees with 306 levels between redshift z — 6 and z = 2. With such a fine time resolution, these 'splitters' become a 
real problem, which needs to be addressed. More specifically, since mergers result in starburst, we have to be sure 
that mergers identified by the FOF algorithm are real. A merger that immediately follows or is immediately followed 
by a splitter might not be real. 

We have designed an algorithm that scans our merger trees and eliminate splitters that appear to be spurious. We 
start by defining a time of stability i s t a b, which is the minimum time a cluster must exist in the simulation to be 
considered a real cluster. When a cluster forms by a merger, and then splits, we calculate the time t sur during which 
the cluster survived before splitting. If t sur < t sta h, we consider that the cluster was not present long enough to be 
considered real, and therefore the merger that formed that cluster and the subsequent splitter that destroyed it were 
both spurious. This handles the first case described above, that is of two clusters that briefly come into contact during 
a non-merging encounter. 

This leaves the second case, two clusters that experience a series of mergers and splitters before finally merging for 
good. In this case, every time the cluster splits, the two pieces will re-merge after a time i rc t- If this time is sufficiently 
short, we assume that the split and the subsequent re-merging were spurious, and that the cluster actually never split. 
We use the criterion £ ret < t s tab, where i s tab is the time of stability that was used in the other criterion. 

We have applied this method to the merger tree built from our simulation, using a time of stability i s tab = 1-5 x 
10 8 years. We have identified all cases where a cluster splits. For each case, we calculated the time t SUT during which 
the cluster existed before it split, and the time t re t elapsed before the various pieces re-merged together (this can be 
infinite if the pieces never re-merged). In Figure [ATI we show a plot of t rct vs. i sur , where the contours show the actual 
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t sur (years) 



Fig. Al. — Contour plot showing the time t Bur that a merger remnant survives before it splits, vs. the time t re t that a split lasts before 
the various parts re-merge together. The contours show the number of clusters. All clusters have either t BUI < t s tabi *ret < t B tabi or both, 
where t sta b = 1.5 X 10 8 years is the time of stability. 



number of clusters. Most clusters have either t sm < t s t a b) ^rct < ^stab, or both. Hence, according to our criterion, these 
splitters are not real, and we simply ignore them. A small fraction of splitters fail both criteria, and are considered real 
splitters. Such splitters occur when a small halo is torn apart by the tidal field of a larger one. This violent process 
should result in a significant reheating of the gas, which can later cool down and lead to a starburst. Consequently, 
we treat a halo formed by a splitter like we treat halos formed by monolithic collapse and mergers: we calculate its 
cooling time, and after that time has elapsed, we start a new outflow. Note that the details of this treatment are not 
critical, since real splitters are rare. 

THE SPHERICAL COLLAPSE MODEL 

In SB01 and Paper I, the spherical top-hat collapse model was used to estimate the size of uncollapsed halos when 
these halos are hit by outflows. The comoving radius of a halo is given by 

/2\ 1/3 JZ(l-CQBg) 

\9J sin 0) 2 / 3 ' 1 ' 

where R is the comoving radius in the limit z — > oc. In equation (|B1[) . the parameter 8 is obtained by solving the 
following transcendental equation, 

-sinfl\ 2/3 _ Mg) 
2vr J 5 + (z con )> 

where S + is the linear growing mode. As z varies from z — oo to z — z co ii, varies from to 2-7T. The physical radius 
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of the halo is given by 



( 



2 



) 



1/3 



R(l - cos 9) 



9 



(i + z)(e-sme) 2 / 3 ' 



(B3) 



Notice that at high redshift (when ft w 1 and A « 0), S + (z) « (1 + z)^ 1 . Hence, the denominator in equation (IB3|) is 
nearly constant. The time-dependence in equation (|B3[) essentially comes from the factor (1 — cos#) in the numerator. 
The halo has a physical size of at 8 = (corresponding to the big bang), reaches maximum size at = tt, and 
collapses to a black hole at 9 = 2ir. 

In practice, halos will not collapse to a black hole. After turnaround, radial motions will be converted into transverse 
motions, and eventually the halo will virialize, with a physical radius that is 1/2 of the maximum physical radius (that 
is, the halo shrinks by a factor of 2 in radius after turnaround). We do not have a precise model for the evolution of 
the system between turnaround and virialization. A fair assumption is that the time it takes for the halo to shrink 
by a factor of 2 is equal to the time it woul d have taken to collapse to a black hole in the absence of virialization 
(|PadmanabhanlfT993l; iDekel fc Ostrikerlll999l see however IColes fc Lucchiniri995| ). To achieve this result, we modify 
equation (|Blj) and (|B3|) . using 



Figure IBTI shows the evolution of a spherical halo collapsing at redshift z co \\ = 3 in a universe with tt = 0.268 and 
A = 0.732. The comoving radius and physical radius are plotted versus time in units of the Hubble time 1/H . Since 
the evolution is similar for halos of different sizes, we set arbitrarily R = 1. The solid curves show the basic spherical 
model given by equations (|B1|) . (|B2[) and (|B3[) . The dashed curves show the effect of using equation (|B4|) and (|B5|) . 
After turnaround, the halo collapses by a factor close to 2. 

It is convenient to re-express the physical radius in terms of its final value at z = z co ih That value is obtained by 
setting z — z co a and 9 = 2ir in equation (IB5[) . We get 




(B4) 



(B5) 



where 




(B6) 



(1 + Z C oll)^+( z coll) 

(l + z)6+(z) 



1 2 /3 




(B7) 



Since at high redshift the quantity in brackets is close to unity, the time-dependence is essentially given by the 
function x(#)- 
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t/H 

Fig. Bl. — Evolution of a spherical halo collapsing at rcdshift z co \\ = 3 in a universe with f2o = 0.268 and Ao = 0.732. The solid and 
dashed curves show the standard and modified spherical collapse model, respectively. 
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